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Abstract 

The focus of the article is on the analysis of a semi-infinite crack at the interface between 
two dissimilar anisotropic elastic materials, loaded by a general asymmetrical system of forces 
acting on the crack faces. Recently derived symmetric and skew-symmetric weight function 
matrices are introduced for both plane strain and antiplane shear cracks, and used together 
with the fundamental reciprocal identity (Bctti formula) in order to formulate the clastic 
fracture problem in terms of singular integral equations relating the applied loading and 
the resulting crack opening. The proposed compact formulation can be used to solve many 
problems in linear elastic fracture mechanics (for example various classic crack problems in 
homogeneous and heterogeneous anisotropic media, as piezoceramics or composite materials). 
This formulation is also fundamental in many multifield theories, where the elastic problem 
is coupled with other concurrent physical phenomena. 

Keywords: Interfacial crack, Stroh formalism, Weight functions, Betty Identity, Singular 
integral. 



1 Introduction and formulation of the problem 

The method of singular integr al equations in linea r elasticity was first developed for solv- 
ing two-dimensional problems, ( Mushkelishvili . 19531 ). and later extende d to three-dimensiona l 



cases by means of multi-di mensional singular integral operators theory (jKupradze et all I197S 



Mikhlin fc Prossdor i flQSnh . Singular integral formulations for both two and three-dimensional 
crack pr oblems have been derived by means of a general approach based on Green's function 
method (jWeaverl . 119771 ; iBudiansky &; Ricd . Il979l : LLinkov et all 119971 ). As a result, the displace- 
ments and the stresses are defined by integral relations involvi ng the Green's functions, for 
which explicit expressions are required (|Bigoni k Capuanil . l2002l h Although Green's functions 
for many two and three- dimensional crac k prob l ems in isotrop i c and anis o tropic elastic mate- 
rials have been derived (jSinclair & Hirthl . 1 19751 : IWeaverl . Il977l : IPanl . l2000l . 120031 : IPan fe Yuan! 
200(1 ). their utilization in evaluating physical displacements and stress fields on the crack faces 
implies, especially in the anisotropic case, challenging numerical estimation of integrals which 
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convergence should be asserted carefully. Moreover, the approach based on Green's function 
method works when the tractions applied on the discontinuity surface are symmetric, but not 
in the case of asymmetric loading acting on the crack faces. 

Recently, using a procedure based on Betti's reciprocal theorem and weight functional an 
alternative method for deriving integral identities relating the applied loading and the resulting 
crack opening has been developed for tw o and three-dimensiona l semi -infinite interfacial cracks 
between dissimilar isotropic materials by iPiccolroaz &; Mishurisl (j2013h . In the two-dimensional 
case, the obtained identities contain Cauchy type singular operators together with algebraic 
terms. The algebraic terms vanish in the case of homogeneous materials. This approach avoids 
the use of the Green's functions without assuming the load to be symmetric. 

The aim of this paper is to derive analogous integral identities for the case of semi-infinite 
interfacial cracks in anisotropic bimaterials subjected to two-dimensional deformations. 

General expressions for symmetric and skew-symmetric weight functions for interfacial cracks 



in two-dimensional anisotropic bimaterials have been recently derived by Morini et al. ( 20131 ) 



by means of Stroh representation of displacements and fields (jStrohl . I1962T) c ombine d with a 
Riemann-Hilbert formulation of the traction problem at the interface ( Suo . 1990bl ). These 
expressio ns for the weight funct i ons ar e used together with the results obtained for isotropic 
media by Piccolroaz Sz Mishuris ( 20131 ) in order to obtain integral formulation for interfacial 
cracks problems in anisotropic bimaterial solids with general asymmetric load applied at the 
crack faces. 

We consider a two-dimensional semi-infinite crack between two dissimilar anisotropic elastic 
materials with asymmetric loading applied to the crack faces, the geometry of the system is shown 
in FigJTJ Further in the text, we will use the superscripts W and ^ to denote the quantities 
related to the upper and the lower elastic half planes, respectively. The crack is situated along 
the negative semi-axis x\ < 0. Both in -plan e and antiplane stress and deformation, which in 
fully anisotropic materials are coupled (jTingj . Il995l ). are taken into account. The symmetrical 
and skew-symmetrical parts of the loading are defined as follows: 



(P) = 2 (P + + P 



1PJJ 



P 



P 



1] 



where p + and p _ denote the loading applied on the upper and lower crack faces, X2 = + and 
X2 = 0~, respectively (see Fig. [IJ). 

In Section [2] preliminary results needed for the derivation of the integral identities and for 
the complete explanation of the proposed method are reported. In Section [2~!T1 the fundamental 
reciprocal identity and the weight functions, defined as special singular solution of the homo- 
geneous traction- free problem are introduced. In Section 12.21 symmetric and skew-symmetric 



Defined by iBuecknerl (|l985l ) as singular non-trivial solutions of the homogene ous traction-free problem and 
later deriv ed fo r general three-dimensi onal problems by IWillis fc Movchanl (|l995h , and for interfacial cracks by 
iGad (|l992h and IPiccolroaz et all |2009l ). 



2 



weight funct i ons m atrices for interfacial cracks in anisotropic bimaterials recently derived by 
Morini et all (|2013h are reported. 



Section [3] contains the main results of the paper: integral identities (|3"3|) . ([35]) . (fM|) and ([65]) 
for two-dimensional crack problems between two dissimilar anisotropic materials are derived 
and discussed in details. The integral identities are derived for monoclinic-type materials, which 
are the most general class of anisotropic m edia where both in-plane and antiplane strain and 
in-plane and antiplane stress are uncoupled fanel Il99fil l2000h , and the Mode III can be treated 
separately by Mode I and II. By means of Betti's formula and weight functions, both antiplane 
and plane strain fracture problems are formulated in terms of singular integral equations relating 
the applied loading and the resulting crack opening. 

In Section [H the obtained integral identities are used for studying cracks in monoclinic 
bimaterials loaded by systems of line forces acting on the crack faces. The proposed examples 
show that using the identities explicit expressions for crack opening and tractions ahead of the 
crack tip can be derived for both antiplane and in-plane problems. These simple illustrative 
cases demonstrate also that the proposed integral formulation is particularly easy to apply and 
can be very useful especially in the analysis of phenomena where the elastic behaviour of the 
material is coupled with other physical effects, as for example hydraulic fracturing, where both 
anisotropy of the geological materials and fluid mo tion must be taken into ac coun t. 

Finally , in A ppendix A, the Stroh formalism ( Stroh . 19621 ). adopted by Su 3 (|l990bh and 
Gao et alJ (1 1992h in a nalysis of interfacial cracks in anisotropic bimaterials and recently used 



by iMorini et al.l ( 20131 ) for deriving symmetric and skew-symmetric weight functions, is briefly 
explained. In particular, explicit expressions for Stroh matrices and surface admittance tensor 
needed in weight functions expressions associated to monoclinic materials are reported. 



2 Preliminary results 

In this section relevant results obtained by several studies regarding interfacial cracks are re- 
ported. These results will be used further in the paper in order to develop an integral formulation 
for the problem of a semi-infinite interfacial crack in anisotropic bimaterials. 

In Section [2TT1 we introduce the Betti integral formula for a crack in an elastic body subjected 
to two-dimensional deformations with general asymmetric loading applied at the faces. 

In Section 12.21 general matrix equations expressing weight fu nctions in terms of the associate 



singular traction vectors, recently derived by lMorini et alJ (120131 ). and valid for interfacial cracks 



in a wide range of two-dimensional anisotropic bimaterials are reported. 
2.1 The Betti formula 

The Betti formula is generally used in linear elasticity in order to relate the physical solution to 
the weight function which is defined as special singular solution to the homogeneous traction-free 
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Figure 1: Two-dimensional semi-infinite interfacial crack loaded by non necessarily symmetric forces applied on 
the crack faces. 



problem (jBuecknerl . Il985l ; IWillis Movchanl . Il995l ). Since the Betti integral theorem is inde- 
pendent of the specific elastic constitutive relations of the material, it applies to both isotropic 
and anisotropic media in the same form. 

The notations u = (u\,U2 : U3) T and rr = (021, 022 > 0"23) T are introduced to indicate respec- 
tively the physical displacements and the traction vector acting on the plane X2 = 0. According 
to the fact that two-dimensional elastic deformations are here considered, both displacements 
and stress do not depend on the variable X3. Nevertheless, since both in-plane and anti -plane 
strain and stress are considered, non- zero components U3 and 023 are accounted for (jTingl . 
19951 ). The notations U = (U u U2, U?,) T and £ = (£21, ^22, ^2^) T are introduced to indicate 



the weight function, defined by iBuecknerl (119851 ) as a non-trivial singular solution of the homo- 
ge neous traction- free p roble m, and the associated traction vector, respectively. As it was shown 
by IWillis &: Movchan (|l995h . the weight function U is defined in a different domain respect to 
physical displacement, where the crack i s placed along the p ositiv e semi-axis x? > . Follo wing 
the procedure reported and discussed in Willis fc Movchan ( 1995 ); Piccolroaz et al. ( 20091 ) and 
Piccolroaz fc Mishurisl ( 2013 ). from the application of the Betti integral formula to the physical 
fields and to weight functions for both the upper and the lower half-planes in Fig. [TJ we obtain: 

r°° ( ~ 

/ {RU(xi -xi,0 + ) -p + (xi) -RU(xi -xi,0~) -p~(xi)+ 
J— 00 

+RU(ii - xi, 0+) • cr(+)(xi,0 + ) - RU(x; - xi, - ) • <r (+) (xi, - )- 
RS(x' 1 -xi,0 + )-u(xi,0 + ) + R£(x / 1 -xi,0~)-u(xi,0") \ dx x = 0, (2) 
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where R is the rotation matrix: 



R 




p + (xi) and p~(xi) are the loading acting on the upper and on the lower crack faces, respectively, 
and cr( + \xi,0) is the physical traction at the interface ahead of the crack tip. The superscript 
( + ) denotes a function whose support is restricted to the positive semi-axis, x\ > 0. In eq. ([2]), 
x'i denotes a shift of the weight function within the plane (xi,X3) and the dot symbol stands 
for the scalar product. 

Assuming perfect contact conditions at the interface, which implies displacement and traction 
continuity at the interface ahead of the crack tip, <j( + ) (xi,0) can be defined as follows: 



r(+) 



(X!,0 H 



r(+) 



-(+) 



(xi), x 1 > 0. 



(3) 



Similarly, also the traction corresponding to the singular solution U satisfies the continuity at 
the interface: 

S(xi,0+) = E(xi,0~) = S(xi), xi < 0. 



(4) 



Using these definitions, (J2J) becomes: 

|R[[U]]K - Xl ) • t«(xi) - RSK - Xl ) ■ M H (zi)} dxi 



f |RU(a;i -xi,0 + ) ■p + (x 1 ) - RU(xi -xi,0 _ ) • p~(a;i)} ctei, 

J— oo ' 



(5) 



where [[u]]( ) is the crack opening behind the tip, ( ' denotes that its support i s restricted to the 
negative semi-axis, x~\ < 0, and [[Uj| is known as t he symmetric weight function (jWillis fc Movchanl . 
19951 : IPiccolroaz etaP . 120091 : iMorini et all l2013h : 



[[U]](x 1 ) = U(x 1 ,0 + )-U(x 1 ,0-). 



(6) 



By expressing the loading acting on the crack faces in terms of the symmetric and skew- 
symmetric parts defined by (H|), the Betti identity ([5]) finally becomes: 



{r[[u]]K 



Xl) ■ T^ + \x\) 



RS 



{R[[U]](x; - xi) • (p) (x!) + R(U) (xl - Xl ) • M&i)} dx! 



(7) 



where (U) is known as the skew-symmetric weight function ([Willis &: Movchanl . ll995l : lPiccolroaz et al. 
20091 : IMorini et aD . \20l$ ) : 



(U) (x 1 ) = -[U(x 1 ,0+)+U(x 1 ,0-)] 



(8) 
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The integral identity (|7|) can be written in an e quiva lent form using the convolution respect to 
x\, denoted by the symbol * (jArfken fc Weberl . 120051 ): 



Capj] 



* T 



(+) 



(RS) T * [[u 



(R[[ur *<p>-(R<u>r *[[ P ]]. 



(9) 



This integral identity relates physical traction and crack opening to weight functions and load 
applied at the crack faces, and will be used further in the text in order to formulate the interfacial 
crack problem between dissimilar anisotropic materials in terms of singular integral equations. 

Note that, in order to simplify notations, in eq. ([9]) the scalar product a • b between vectors 
a and b is replaced by the "row by column" product a r b between the row vector a T and the 
column vector b. 

2.2 Symmetric and skew-symmetric weight functions for anisotropic bimate- 
rials 

Let us introduce the Fourier transform of a generic function / with respect to the variable x\ 
as follows: 

/oo 1 /•oo 

f( Xl )e^d Xl , f( Xl ) = ^ [/(£)] = - / f(Oe-«*idZ. (10) 
-oo £K J — oo 



In iMorini etaD (120131 ). the following expressions for the Fourier transform of the singular dis- 
placements U at the interface between two dissimilar anisotropic media have been derived: 



u(e,o + ) = {^(yw-y (1) ) - ^(y« + y«)}s-(o, e € R, 

U(£,(T) = {1(Y( 2 ) - Y (2) ) + ^(Y( 2 ) + Y (2) )}E"(0, e € R, 



fill 



where 5] is the Fourier transform of the singular traction at the interface, which in the case 
of perfect contact condition is de fined as in expre ssion and Y is the Hermitian definite 
positive surface admittance tensor ()Gao et all Il992l ). depending on the elastic properties of the 
materials and defined in details in Appendix A. 

The superscripts + and ~, used here and in the sequel, denote functions analytic in the upper 
and in the lower complex half-planes, respectively 



Eqs. (|lip and (I12j) represent general expressions relating the tractions 5] (£) applied on the 
bounding surfaces and the corresponding displacements U(£,0 + ), U(£, 0~) for the upper and 
lower half-planes, respectively. 
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The symmetric and skew-symmetri c weight function matrices are derived by taking respec- 
tively the jump and the average of TJ ([Willis &: Movchanl . Il995l ; iPiccolroaz et al.l . l2009l ): 



Eqs. (USD and £E 



[[U]] + (0 = sign(0 Im(Y( 1 ) - Y( 2 )) - Re(Y« + Y^)}t (£), 

^{i sign(e) Im(Y« + Y ( 2 )) - Re(Y« - Y( 2 ))}sT(£), £ £ R. 
can also be expressed in the compact form: 

[[U]] + (0 = "T^{ ReH " * si ^(0 ImH}s~(e), 



L{ ReW -isign(£) ImWjs (0, (GR, 



2iei 



where H and W are the bimaterial matrices defined as follows JSu3, Il990bl : iTinel l2000h : 



H 

W 



yW+Y 

yW-y 



(2) 
(2) 



(13) 
(14) 

(15) 
(16) 



(17) 
(18) 



Expressions (|15p and f)16|) are valid for interfacial cracks in general anisotropic two-dimensional 
media. Since in anisotropic materi als in-plane and antiplane displacements and stresses are 
generally coupled jXineL Il99d . l2OO0h . for the case of fully anisotropic media three linearly inde- 
pendent vectors S and then U must be defined for obtaining a complete basis of the singular 
solutions space. Nevertheless, there are several classes of anisotrop i c mat erials where in-plane 
and antiplane displacements and stress are uncoupled fanel Il995l . l200d ) and then Mode III 
deformation can be treated separately from Mode I and II as for the case of isotropic media 
( Piccolroaz et al. . 20091 . l20ld : IPiccolroaz fc Mishurisl . \2QlA ) . . n the next Section integral identi- 
ties are derived for interfacial crack between two media belonging to the most general of these 
classes, known as monoclinic materials. 



3 Integral identities 



In this Section, following the approach of Piccolroaz Sz Mishuris ( 20131 ) . an integral formulation 
of the problem of a semi-infinite two-dimensional interfacial crack in anisotropic bimaterials is 
obtained. A particular class of anisotropic materials, where elastic properties are symmetrical 
with respect to a plane, is considered. These materials are known as monoclinic, and in the case 
in which the plane of symmetry coincides with x% = both in-pla ne and antiplane displacements 



1995). Monoclinic having plane of 



and in-plane and antiplane stress are uncoupled ([Ting], 
symmetry at X3 = are the most general class of anisotropic materials where stress and strain 
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are decoupled, and they include as subgroups all other classes hay ing this property, such as 
orthotropic and cubic materials (jHorgan fc Milleri . Il994l : iTinal . \200d ). Explicit expressions for 
Stroh matrices and surface admittance tensor corresponding to these type of media are reported 
in Appendix A. These expressions have been used for evaluating bimaterial matrices (|17[) and 

C5D. 

In Sections 13.11 and 13.21 antiplane shear and plane strain interfacial cracks problems in 
monoclinic bimaterials are formulated in terms of singular integral equations by means of weight 
function expressions (|15p and (|16p and Betti integral identity Q. 



3.1 Mode III 

Considering antiplane deformations in monoclinic materials, as it is shown in Appendix A, 
constitutive relations reduce to scalar equations relating stresses a 23 and a\ 3 to U3, and then the 
traction o 23 = rj~ an d the displacements derivative for both upper and lower half-plane material 
become |Suol . Il990bh : 

(19) 



T3(x 1 ,x 2 ) = L 33 h 3 (z 3 ) + L 33 h 3 (z 3 ), 



u 3 ,i(x 1 ,x 2 ) = F 33 h 3 (z 3 ) + F 33 h 3 (z 3 ), 
where z 3 = xi + ^3X2. The bimaterial matrices (TT71) and (fl~8l) reduce to: 



#33 
W 33 



s 44 s 55 



' 15 



(1) 



+ 



s 44 s 55 



'15 



s 44 s 55 



' 15 



(1) 



s 44 s 55 



'15 



(2) 
(2) 



(20) 

(21) 
(22) 



where Sj- are elements of the reduced elastic compliance matrix (see Appendix A). According 
to general expressions (|15p and (|16p . the Fourier transform of symmetric and skew-symmetric 
weight functions for an antiplane shear crack between two dissimilar monoclinic materials are: 



H- 



33 



^23(0, (U3K0 



~t^) = ^[L> 3 ]] + (0, 



(23) 



where the following non-dimensional parameter has been introduced: 



s 44 s 55 



^15 



(1) 



s 44 s 55 



5 45 



(2) 



• s 44 s 55 



' 15 



(1) 



+ 



s 44 s 55 



'15 



(2)- 



(24) 



In the case of antiplane deformations, the Betti formula reduces to the scalar equation: 

[[u 3 ]] * r« - s 23 * [[« 3 ]](-) = -M * (p 3 ) - (Us) * [[Psl 



(25) 
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Applying the Fourier transform with respect to x±, defined by relation (jlOp . to this identity, we 
obtain: 

M + f 3 + " S 2 - 3 [[n 3 ]]- = -[[^3]] + (pa) - (^s)[[p 3 ]]. (26) 
Multiplying both sides of (|26p by [[?7 3 ]] _1 , we obtain: 

f 3 +-£[Nr = -(P3}-4[P3]. (27) 
The factors in front of the unknown functions are given by: 



A = au 3 ]r i (u 3 ) 



2' 



I^3]]- 1 S 3 



33 



(28) 



If we apply the inverse Fourier transform to (|27p , we derive two distinct relationships corre- 
sponding to the two cases x\ < and x\ > 0: 



(p 3 )(xi) + -F- 1 1 <0 U[N] 



J" 



-l 



xi<0 



x\ < 0, 



t+(x 1 )=T- 1 >0 \bIu 3 ]}-1 xi >0 



(29) 
(30) 



It is important to note that the term fg cancels from (|29p because it is a "+" function, while 
\p 3 ] and (]5 3 ) cancel from (|30j) because they are "— " functions. 

To proceed further, we need to eval uate the inverse Fourier t ransfo rm of the function 
Following the procedure illustrated by Piccolroaz Sz Mishurisl ( 2013 ). we get: 



1 d[[u 3 f 

— * 



irx\ 



dx\ 



1 d[[u 3 ]Y 
x\ — r] drj 



■drj. 



(31) 



Then we can define the singular operator S and the orthogonal projectors V± (P+ + V— = 1) 
acting on the real axis: 



1 1 f°° 

ib = Sip = * ^(^l) = — — 

TTXl 7T J-oo X\ 



dr), 



±xi > 0, 
otherwise. 



(32) 



(33) 



The operator S is a singular operator of Cauchy type, and it transforms any funct ion tp satisfying 
the H older condition into a new function S<p which also satisfies this condition (jMushkelishvilil . 
T he properties of the operator S in several functional spaces have been described in 



details in iProssdora ()1974l ) . 

The integral identities (|29p and (|30p for a Mode III interfacial crack between two dissimilar 
monoclinic materials become: 



(P3>0=i) + jMCzi) 



i 



H- 



d[[n 3 ]] 



(-) 



33 



dx\ 



xi < 0, 



(34) 
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.(+) 



(*l) 



xi > 0, 



(35) 



#33 dx\ 

where = VSV- is a singular integral operator, and = V+SV is a compact integral 
operator ( Gakhov Sz Cherski , 1978 ; Krein, 1958 ; Gohberg fe Krein . Il958l ), These two operators 
look similar, but they are essentially different, in fact: : -F(R_) — > -F(R_), while : 

-F(R_) —> -F(R+), where .F(R±) is some functional space of functions defined on H±. 

For explaining better this point, the integral identities (|34p and ()35[) can be written in the 
extended form: 



i 



vri?33 



1 dlu 3 }}^ 
x\ — 7] dr] 



dr], x\ < 0, 



.(+) 



(*i) 



1 



1 « H 



dr/, x\ > 0. 



(36) 



(37) 



ttH 33 J-oc x\ -7] dr] 

The integral in ()36[) is a Cauchy-type singular int egral with a moving sing ularity, whereas the 
integral in (I37j) possesses a fixed point sing ularity dPuduchaval . Il97fil ll979T ). 

In the case of a homogeneous monoclinic material, the integral identities ()34|) and (|35p 

simplify, since = 0, #33 = 2 J S44S55 — S45 and thus there is no influence of the skew-symmetric 
loading. 

Summarizing, the integral identities for Mode III interfacial cracks in monoclinic bimaterials 
are given by equations (|34[) and (|35p . The equation (|34|) in an invertible singular integral 
relation between the applied loading (^3), [[p 3 ]] and the corresponding crack opening [[143]] (~\ 
The equation (|35p is an additional relation through which it is possible to define the behaviour 
of the solution [["Us]]^ - -*. Since the operator S^ c ' is compact, it is not invertible, and thus for 

tip 

the equation ([34"|) (see iMushkelishvilil (|l946l ) for details). 



deriving the traction ahead of the cr a ck tip r^f' one needs to evaluate [[1*3]]^ ^ by inversion of 



3.2 Mode I and II 

For plane strain deformations i n monoclini c materials, the surface admittance tensor Y is given 
by a 2 x 2 matrix of the form ( Ting . 19951 ): 



where: 



Y = s u P + i(s n c- s 12 )E 
'0 -1 



b d 

d e 



E 



1 

Ml + M2 = a + ib, 11x^1 = c + id, 
e = ad-bc = Im[/ii^ 2 (Pi 



(38) 

(39) 

(40) 
(41) 
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in which fj,\ and fi2 are solutions of the eigenvalue problem as sociated to balanc e equa tions 
by means of Stroh representation of displacements and stresses ( Stroh . 19621 ; Tine . 19961 ) (see 
Appendix A for more details). Since \i\ and [12 are eigenvalues with positive imaginary part, b 
is strictly positive, b = Im(/ii + H2) > 0, w hile the positive definiteness of the matrix Y, and 
consequently of P, implies that (jTingi . Il995l ): 



e > and be - d 2 < 0. 



(42) 



Thus, bimaterial matrices H and W for an interfacial crack between dissimilar monoclinic 
materials under plane strain deformations can be decomposed into real and imaginary parts as 
follows: 



"(1) 



+ Y 



(2) 



H + i^H u H 22 B, 



H = 

W = Y (1) - Y {2) = W' - Hv/tfii^E, 
where matrices H' and W' are defined as: 



H 



«V H11H22 

H22 



w 



\ ^VH\lH22 



\\[RxyH^2 
S2H22 



(43) 
(44) 

(45) 



Note th at Hu and H22 are real positive parameters defined similarly to those introduced by lSuo 
dl990bl ) for orthotropic bimaterials: 



H11 = + KiP, H 22 = [es n ]^ + [es 



1(2) 



(1) 



11J 



1(2) 



(46) 



Regarding matrix H, two non -dimensional Dundurs-like parameters are defined (ITingl . Il995l : 
SuflL Il990bl : iMorini et all l2013h : 



a 



[cfaljW + [d S ' n ] (2) , p = l£n£_- 4] (1) " [4c- S ' 12 p) 



\JHwH22 



\fH\xH2f2 



(47) 



while the matrix W depends by four non-dimensional Dundurs-like parameters (ITingl . Il995l : 
Suol . Il990bl : IMorini et all 120131 ) : 



A 



[bs' n ]W - [bs' n ]W 

61 = 

[ds' u ]M - [efaiJP) 



r.s 



11 



(1) 



r.s 



11 



(2) 



\fHx\H22 \/H\\H22 

The Fourier transforms of the symmetric and skew-symmetric weight functions (|15p and (JX6 
for a plane monoclinic bimaterial assume the form: 



H22 

[siiC-s' 12 ]W + [4c-4]^ 



(48) 



(49) 



[[U]] + (0 = - i (H - isign(e)/3v^iT^E) ± (0, 



(50) 
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(U)(e) = -^(w'+i S ign(e) 7 V^^E)s (£). (51) 

Since in plane strain elastic bimaterials Mode I and Mode II are coupled, two linearly indepen- 
dent singular solutions and tractions LP = (Ul^U^) 7 ', S* = (Y;\-i ; Ti\ ) T ,i = 1,2, are ne eded 
in order to define a complete basis of the singular solutions space (jPiccolroaz et all 120091 ). As 
a consequence, in this case symmetric and skew-symmetric weight functions [U] and (U), and 
the associate traction I] are represented by 2 x 2 tensors which may be constructed by ordering 
the components of each singular solution in columns: 

U=(3 E*V 90. (52) 



U 2 U 2 J V ^22 ^22 



Correspondingly, the rotation matrix R reduces to: 



;). (53) 



Applying the Fourier transform to the ([9]), we obtain: 

[[U]f Rr+ - S T R[[u]]- = -[[UfR(p) - <U) T R[[p]], £ £ R. (54) 
Multiplying both sides by R -1 [[U]] _T , the following identity is derived: 

f + -B[[u]]- = -(p)-A[[p]], (55) 

where A and B are given by: 

A = R" 1 [[U]]- T (U T )R, B = R- 1 [[U]]- T S T R. (56) 



Explicit expressions for these matrices can be computed using symmetric and skew-symmetric 
weight functions (f50l) and (f5"Tj) : 

(A'+isign(£)A"), (57) 



2^H n H 22 (a 2 + /3 2 -l) 



B = 7W tt I^mR2 n ( B ' + i^COE), (58) 
^H u H 22 (a 2 + fi z - 1) 



where A , A and B are: 



A' = f VHnH 2 2(a\ - /?7 - 5i) ff 22 (A-a<5 2 ) 

ffuCA-aJO ^H u H 22 (aX- Pj - 6 2 ) 



^H u H 22 ( ai + (3X) g 22 ( 7 + /3(5 2 ) 
-#n(7 + /%) V^n^2 2 (a7 + /3A) 



(59) 
(60) 
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Applying the inverse Fourier transform to the identity (|55p . for the two cases x\ < and x\ > 
we get: 

x\ < 0, 



(p) (xi)+ J Xl 1 <0 [A[p]_ 
r(+)(^ 1 )+J-- 1 >0 [A[p] 



• r x 1 <0 
- • r x 1 >Q 



B[u]" 
B[ul 



xi > 0. 



(62) 
(63) 



As for the case of antiplane deformations, illustrated in the previous Section, the term r + in eq. 
(]55p cancels from the (I62D because it is a "+" function, while (p) cancels from the (I63p because 
it is a "— " function. 

Using the same inversion procedure of the previous Section the following integral identities 
for plane strain deformations in monoclinic bimaterials are derived: 



(p) (x 1 ) + A^M=^ ( ' ) — 



-) 



\x 1 ) + A^[[p]]=B 



(c) 



d[[n]](~) 



dxi 



xi < 0, 



x\ > 0, 



(64) 



(65) 



where matrix operators A {s) ,B^ s) : F(R_) -»• F(R_), and A {c \B {c) : F(R_) -»• F(R + ) are 
defined as follows: 

-(A' + A"S (S) ), (66) 



2VH u H 22 (a 2 + /3 2 - 1 
1 



V#ntf 2 2(a 2 + /3 2 -l) 
1 



2^H u H 22 (a 2 + /3 2 - 1 



(b'S^ -/3E), 
A"5W, 



S (e) 



(67) 
(68) 
(69) 



^H u H 22 (a 2 + /3 2 - 1) 

Equations ([64^1 ([65]) . together with the definition of operators ([66]) . ([67]) . ([68]) and ([69]) . form 
the system of integral identities for Mode I and II deformations in monoclinic bimaterials. The 
equation (|64p is a system of two coupled singular integral equations, which decouples in the 
case where the Dundurs parameters a and (3 vanish. Observing expression (|47p . we can note 
that (3 vanishes in the case of a homogeneous monoclinic material, while a is zero only for 
some particular subclasses of materials, such as for orthotropic materials, where the quantity d, 
defined by (I40p and representing the imaginary part of the product of the eigenvalues, vanishes 
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jSixoL Il990al : iGupta et all 119921 1. As a consequence, for a homogeneous orthotropic material, 
the system (j64"|) is reduced to the following decoupled equations: 



H 



1 S {s)dln 1 f-) 
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H- 



1 s ts)dh 2 f-) 



22 



(pi) Oi 

(P2 




S^M, xi<0, 



t/^WM, xi <o. 



(70) 
(71) 



The solution of these two equations requires the inversion of the singul a r ope rator S^ s \ which 
has been performed and discussed in details by Piccolroaz fc Mishuris ( 20131 ). The inversion 
of the matrix operator in the general case requires the analysis of the systems of singular 
integral differential equations ( Vekua . 1970l ). 



4 Illustrative examples: line forces applied at the crack faces 

In this section we report an illustrative example of application of the integral identities in analysis 
of interfacial cracks in anisotropic bimaterials. Antiplane (Fig. [2| and plane strain (Fig. [3|) 
interfacial cracks in monoclinic bimaterials loaded by line forces acting on the crack faces are 
studied by means of the proposed integral formulation. Explicit expressions for crack opening 
and tractions ahead of the tip corresponding to both symmetrical and skew-symmetrical loading 
configurations are derived. The proposed illustrative cases show that the integral identities 
derived in previous Section represent a very useful tool for studying interfacial crack problems 
in anisotropic materials, and their relevance is even greater in analysis of scenarios where the 
elastic problem is coupled with other phenomena, as for example hydraulic fracturing where 
both anisotropic elastic behaviour of geomaterials and fluid motion must be taken into account. 
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(A) Mode III: symmetrical forces (B) Mode III: skew-symmetrical forces 



Figure 2: Geometry of the Mode III crack. 
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(A) Mode I and II: symmetrical forces 



(B) Mode I and II: skew-symmetrical forces 



Figure 3: Geometry of the Mode I and II crack. 

4.1 Mode III: symmetrical line forces 

We consider an antiplane shear crack where the loading is given by two symmetrical line forces 
applied on the faces, at a distance a from the crack tip, and directed along the X3-axis: 



(P3> (xi) = -F6( X1 + a), IbaKsi) = 0, 



(72) 



where 5 is the Dirac delta function, and F has the dimensions of a force divided by a length. 

For antiplane deformations, the singular integral equations relating the applied loading and 
the resulting crack opening is given by (1341) . Inverting the operator by means of procedure 
proposed by iPiccolroaz &: Mishurisl (J2013J), we obtain: 



dim 



l(-) 



H 



7T 



F 



rj 5(t] + a) 
oo V x\ xi-rj 



drj 



H 



7T 



F 



1 



x\ xi + a 



(73) 



Assuming the condition that the crack opening vanishes at zero and at infinity, and integrating 
the ([73]) . the following expressions are derived: 



2F I xi 

[MK^i) = — #33 arctanhW , 

7r y a 



lu 3 ]](xi) = —H33 arctanh*/ — — , 

7T \ Xl 



—a<xi< 0, 
xi < —a. 



(74) 



Substituting ()73|) into equation (|37p . the explicit expression for the traction ahead of the crack 
tip can be evaluated: 



1 



1 gMT 

x\ — rj drj 



■dr\ 



1 



7T V X\ xi + a 



(75) 



It is important to note that the derived expressions for the traction ahead of the crack tip, 
and the associated crack opening, are identical to the results derived by IPiccolroaz fc Mishuris 
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( 20131 ) for the isotropic case, except for the parameter H33, that is related to the anisotropy of 
the media. 

The traction expression (|75p can then be used for calculating the stress intensity factor. 
Applying to this specific case the g e neral d efinition introduced for interfacial cracks in anisotropic 
materials by[Wu| dl99Ch and|Hwuj jl993aJ ) we obtain: 



Km = lim„ \/2vrxiT3 + ' ) (3;i) 



ii->0 



—F. 

7ra 



(76) 



The st ress intensity factor ([76]) , calculated using ([75]) , is identical to that obtained by iHwu 
(|l993al ) for the same symmetric loading configuration applying a different procedure. This is an 
important proof for the validity of the traction expression (|75p and consequently of the associated 
crack opening (|74j> . derived by means of integral identities. 



4.2 Mode III: skew-symmetrical line forces 

We now consider an antiplane shear crack where the loading is given by two skew-symmetrical 
line forces applied on the faces, at a distance a behind the crack tip, and directed along the 
X3-axis: 

(P3> (si) = 0, baK^i) = -2FS(X! + a). (77) 
Applying the inverse operator to equation ([34]) . we obtain: 



H33 j-, 
-v b 



77 5(tj + a] 



drj 



H33 j? 
-v b 



dx\ it J~oo V xi x\ — 77 

Integrating this expression the crack opening is then derived: 

2F 



~a_ 1_ 

x\ x\ + a 



(78) 



X\) = v ^33 arctanh 

7T 

2F 



xi 
a 



-a < x\ < 0, 



[MK^i) = v — #33 arctanh. / , x\ < 

TT V X\ 



Substituting (|T8j) into equation (|37p . the traction ahead of the crack tip becomes: 



+ 



(79) 



(80) 



Th e derived expressions (|79p . (|80p are consistent with the results obtained in lPiccolroaz Mishuris 
(|2013l ^ br isotropic media, the only difference consists in the parameter v which characterizes 
antiplane deformations in anisotropic bimaterials. 

Also in this case, since traction expression (|80l) have been derived, the stress intensity factor 
can be evaluated using the general definition: 



Km 



lim y/2irxiTn + \xi) 



F. 



(81) 
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The explicit expressions for the crack opening (|74p and ()79|) . and for the tractions ahead of 
the tip (|75p and (|80p have been derived only by inversion of the operator and by a simple 
integration procedure. This example shows that the integral identities obtained in Section [3] 
are fairly easy to use for solving antiplane interface crack problems in anisotropic materials. 
Then the integral identities represent a powerful tool for solving many problems in linear elastic 
fracture mechanics, especially for studying systems where the elastic fields are coupled with 
other physical phenomena. 



4.3 Mode I and II: symmetrical line forces 



The plane strain crack problem is now addressed. We assume that the loading is given by two 
symmetrical line forces applied on the faces at a distance a from the crack tip and directed 
respectively along x\ and x 2 axis: 

(pi) (xi) = -FtS( Xl + a), MOcx) = 0, (82) 
<P2> (xi) = -F 2 S(x! + a), [[p 2 ]](xi) = 0, (83) 

where F\ and F 2 have the dimensions of a force divided by a length. 

In two-dimensional interfacial elastic crack problems (plane strain or plane stress), Mode I 
and I I are in general coup led and stress oscillations near the crack tip must be taken into account 
{5ug1 . 1 1990b! : Illwul . Il993al ) . The corresponding integral identities (|64j) are a system of equations 
coupled by means of the generalized Dundurs par ameter f3 (see equation (|4T|) ). connected to the 
oscillation index of the bimaterial e |Suol . ll99nbh : 



2vr 



In 



1-/3 
1 + P 



(84) 



For simplicity, we assume that f) is zero, this implies that also the oscillation index (184ft vanishes 
and that in this particular case we have no oscillations at the crack tip. The equations (I64j) for 
X\ < become: 

^ 5(s) d[M« 



(pi) 

(P2) 



(xi) 



1 

VH n H 22 (a 2 -T) 
1 

VH n H 22 (a 2 -T) 



H 



ii 



a 



d[[ ul 



dx\ 

(-) 



+ 



dx\ 



+ 



S (s) 



22 



gM 

dxi 
dxi 



(-) 



(-) 



(85) 
(86) 



Applying the inverse operator (S^) 1 to both these equations, by means of some algebraic 
manipulations the following result is obtained: 



3[M H . 




dxi 


7T 




H 22 


9xi 


7T 




(87) 



x\X\+ a 
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Then, after integration, for —a < x\ < we derive: 



\u x 



\(-\ Xl ) = 2J ^Ul- Wf^l arctanh^ , 



7r \ V H 2 2 / V a 



[^2]]^ ^(^l) = — \ F2 — aJ—^-Fi ] arctanfu 



And for xi < —a: 



[M (_) (zi) = ^ Ufe - «y f^^ij arctanh^J. (90) 
The tractions components ahead of the crack tip can be evaluated by (1651) : 



^> ^ <«> 

We can observe that, also for the two-dimensional vector problem, if the loading is given by 
symmetric line forces applied at the faces, the corresponding expressions for the traction com- 
ponents ahead of the crack t ip and for stress intensity factors are analogous to that derived by 
Piccolroaz fc Mishurid (|2013h for isotropic bimaterials, the only difference consists in the con- 



stants. The crack opening components (|89j) and (|90j) are coupled to F\ and F2 by means of 
the Dundurs parameter a. This means that, even in the case where the symmetric load pos- 
sesses only one non-zero component, directed along xi-axis or X2-axis, the induced crack opening 
has both x\ and X2 components. This aspect is an important difference respect to the case of 
isotropic materials, and it is connected with anisotropic properties of monoclinic m edia . 



In the particul ar case where a = 0, corresponding to orthotropic bimaterials, (|Sud . Il990al : 



Gupta et all 119921 1. crack opening expressions are similar to the isotropic case, except for the 



different constants. For —a < x\ < 0: 



\ui\y- 'iii = Fiarctanhi/ 

7r V a 



And for x\ < —a: 



I«2]] H (zi) = — F 2 arctanh,p^. (93) 
7r V a 



2H n „ . , I — cT 



\ui]]( Hxi) = Fiarctanh. 

7T V X\ 
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[U2 



l(-) 



2-^22 



i^arctanh. / . 

7T V xi 



(94) 



The tractions expressions (|91 j) and (|92p can then be used for evaluatin g th e stress inten sity 
factors. Applying the general definition explained in details in IWu (jl990h and[Hwu| (|l993al l to 
this particular two-dimensional case without oscillation, we obtain the vectorial formula: 



K = lim V2nx^T {+ > (xi) 
x\ — 



(95) 



where K = (Kjj, Kj) and 
limit, we get: 



-(+) 



(+) J+hr 



Inserting (|9ip and (|92p and evaluating the 



F 2 , K n 



(96) 



Also in this case, as is it has been detect ed for the an tiplane problem, the derived stress 
intensity factors are equal to that obtained by iHwul (|1993al ) in the limit of vanishing oscillatory 
index. This demonstrates that the expressions (|91|) and (|92|) . obtained by means of the integral 
identities, are correct. 

In several studies, a unique complex stress intensity factor accounting both Mode I and 
Mode II c ontributions is associated to two-dimensional int erfacial crack problems in anisotropic 
materials ( Suo . 1990bl : Suo et al. . 19921 : Morini et al. . 2013 ) , similarly to what happens due to the 



oscillations for cracks betwe en different isotropic me dia (IMantic fc Paris! . 120041 : iPiccolroaz et al 
20071 : iGraciani et all 120071 ). It has been shown bv IHwul (I1993ar ) that complex stress intensity 



factor, introduced for anisotropic bimaterials by ISuo ll990bl ). is linked to the classical stress 
intensity factors by a vectorial transformation. In our example, we have assumed that the 
generalized Dundurs parameter (3 and the correlated oscillatory index are zero. In this partic ular 
case, w here there a re no o scillations at the interface and Mode I and II can be decoupled 



Suo 



(|l990bh and IHwul (jl993ah analyses lead to the same expressions (|95|) for real stress intensity 
factors. 



4.4 Mode I and II: skew-symmetrical line forces 

We assume that the loading is given by two skew-symmetrical line forces applied on the faces 
at a distance a from the crack tip and directed respectively along xi and x 2 axis: 

(pi) (xi) = 0, [bi]](xx) = -2FiS(xi + a), (97) 
(pa) (xi) = 0, Mfci) = -2F 2 8(xi + a). (98) 

Also in this case we consider the case where the Dundurs parameter j3 associated to crack tip 
oscillations vanishes. Furthermore, also A is assumed to be zero. Applying the inverse operator 
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(S^)" 1 to the (J63D for x\ < we obtain: 

d[[ Ul ]]^ f, u r-ir 1 



<SifW ■ 1 F 2 VH n H 22 S(x 1 + a), (99) 



3xi 71" V ^1 + a 



S 2 H 22 , + jF^HnH^ix! + a). (100) 



OTl 7T V Xl + a 

Integrating these expressions, for —a < x\ < we derive: 



And for xi < —a: 



ktlF (xi) = <5i-ffiiarctanhw , 

7r V a 



[ U2 ]](-)(xi) = ^k«5 2 if 22 arctanh,/^i, (101) 
7r V a 



2Fi 



[kti]] 1 - Vxi) = — i-^iffnarctanh./ — — + r yF 2 \/H n H 22 , 

7T \ Xl 

= — £ 2 #2 2 arctanh</^- - 7 FWH U H 22 , (102) 

7T V ^1 

Ex cept for the constants, t hese e xpressions for the crack opening are identical to those obtained 
by Piccolroaz Sz Mishurisl ( 2013 ) for the isotropic materials and benchmarked by a comparison 



with the Flam ant solution for a half-plane loaded by two concentrated forces at its surface 
(jBarbeil I2002] ). The tractions components ahead of the crack tip become: 





^2 (+) (^i) = n 1 2 , ( S 2 F 2 + a8 1 F 1 J^ ) J-^—- (104) 
7r(l - a z ) y y H 22 J V xi xi + a 

From these expressions we can observe that, if the loading is given by skew-symmetric line forces 
applied at the faces, the corresponding traction components (|103p and (|104p are coupled to Fi 
and F 2 by means of the Dundurs parameter a. As a consequence, similarly to what we have 
detected for the jump behind the tip in presence of symmetric loading, even in the case where the 
skew- symmetric load possesses only one non-zero component, directed along xi-axis or x 2 -axis, 
the associated tractions have both xi and x 2 components. This coupling effect is not found in 
isotropic materials, and it is connected with anisotropic properties of monoclinic media. 

The stress intensity factors can then be evaluated substituting the tractions ()103|) and (|104p 
into the general expression ([95]): 




Kr = t^xI— \5 2 F 2 + a5 1 F 1 J^ ) , (105) 



Tra \ V ri' 
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K u = T-^\l^- + <^ F 2\ lp-) ■ ( 106 ) 

1 — a z V ira V y iiu y 

For a = 0, the materials assume orthotropic behaviour, and we recover tractions expressions 
similar to those obtained for isotropic media: 

Ti (+) (*i) = %/f (107) 
7T V %i x\ + a 

r< +) (*i) = (108) 
The stress intensity factors become: 

= 5 2 \I-^F 2 , K n = hJ^-F\- (109) 
V 7ra V 7ra 

These expression s for iff and Kjj, corr e spond ing to orthotropic bimaterial, possess the same 

form obtained by Piccolroaz fc Mishuris ( 20131 ) for isotropic media. 



The skew-symmetric stress intensity factors (| 105j) and (|106p have been normalized, multiplying 
respectively by ^/TTa/(^/2F 2 5 2 ) and y/iva/(y/2Fi5i), and plotted in Fig. [Has functions of the 
Dundurs-like parameter a. In order to satisfy the positive definiteness of the bimaterial matrix 
Y, and consequently to verify relations (j4~2j) ( Ting . 19951 ). the value of a must be included in 



an interval a m [ n < a < amax, where the bounds depend on the parameters Hn,H 22 ,5i and 
8 2 . Assuming H\\ = 2.01, H 22 = 6.98, 5\ = 0.72, 5 2 = 0.92, which corres pond s appro ximately 
to a system where material (1) is boron and material (2) is aluminium (Suo, 1990a! ) . we get 



—0.9839 < a < 0.9839. The normalized stress intensity factors computed for values of a within 
this interval and for different values of the ratio F 2 jF\ between the applied forces are reported in 
Fig. 01 As it can be deduced by normalizing expressions (I109P , for a = the curves collapse and 
the materials assume orthotropic behaviour. On the other hand, as the value of a approaches 
the bounds amax = 0.9839 and a m i n = —0.9839, the absolute value of the skew-symmetric 
stress intensity factors increases and consequently the response of the system to skew-symmetric 
loading is amplified. 

The examples illustrated in this Section and in Section 14.31 show that the integral identi- 
ties can be profitably used for solving vectorial interfacial crack problems in two-dimensional 
anisotropic elastic media. The crack opening and the tractions ahead of the tip can be derived 
without the use of Green's function, which implies challenging calculations. For this reason, also 
in the vectorial case, the proposed singular integral formulation may represent a very suitable 
technique in analysis of several fracture processes, especially where coupled fields are involved. 

5 Conclusions 

The problem of a two-dimensional semi-infinite interfacial crack in anisotropic bimaterials has 
been formulated in terms of singular integral equations by means of weight functions and inte- 
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Figure 4: Normalized stress intensity factors corresponding to skew-symmetric forces (see Fig. [3)3) as a function 
of the bimaterial parameter a for —0.9839 < a < 0.9839, computed for Hn — 2.01, H22 = 6.98, Si — 0.72,52 = 
0.92 and different values of the ratio F 2 /Fr. F2/F1 = 0.2, F2/F1 = 0.5, F2/F1 = 1, F 2 /F 1 = 2. 



gral transforms. The proposed method avoids the use of Green's function and the challenging 
numerical calculations related to such approach. Integral identities relating the applied loading 
and the corresponding crack opening have been obtained for both in-plane and antiplane strain 
problems in anisotropic bimaterials. Detailed explicit derivation of the identities have been 
performed for monoclinic materials, which are the most general class of anisotropic media al- 
l owing decoupling between both in-plane and antiplane stress and in-plane and antiplane strain 
f|Tingl . ll995l . l20nd ). As a consequence, the obtained integral formulation can be directly used for 
studying cracks propagation in all anisotropic medi a where in-p lane and antiplane problems are 
decoupled, such as cubic and orthotropic materials (jSuol . 1990b ). Thanks to the great versatility 
of the Stroh formalism, the developed approach can also be easily adapted for studying fracture 
phenomena in many other materials, such as piezoelectrics, poroelastics, and quasicrystals. 

An example of application of the integral identities to antiplane and plane strain crack 
problems where the loading is given by line forces acting on the faces has been performed. 
Explicit expressions for crack opening and tractions ahead of the tip corresponding to both 
symmetric and skew-symmetric loading configurations have been obtained. 

The derived integral identities have various relevant applications especially to multifield 
problems, where the elasticity equations are coupled with other concurrent phenomena, for 
example, but not only, to hydraulic fracturing modelling. Furthermore, they also have their own 
value from the mathematical point of view, as, to the authors best knowledge, such identities 
written in a similar explicit form for interfacial cracks in anisotropic bimaterials seems to be 
unknown in the literature. 
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A Anisotropic materials: Stroh formalism 



In this Appendix, complex variable representation for the stress field and displacements in 
anisotropic e lastic materials subject to two-dimensional deformations based on Stroh formalism 
(IStrohl . ll962h is summarized. Explicit expres sions f or Str oh matrices corresponding to monoclinic 
materials with symmetry plane at X3 = ( Ting . 19951 ) are reported. The surface admittance 
tensor needed in the evaluation of bimaterial matrices (|15p and (|16p and then of weight functions 
([TTD and (HU) is derived. 



A.l Complex variable representation of stress and displacements 

For two-dimensional problems in anisotropic elastic materials, displacements and stress fields 
can be represented in terms of complex variab l e func tions matrices by means of two alternative 
formulations, proposed respectively by Stroh ( 19621 ) and Lekhnitslril ( 19631 ). Introducing the 
stress vectors = (a\k, o"2fc, (?3k) T , k = 1,2 and the displacements u = (u±, U2, u^) T , the consti- 
tutive relations connecting the stresses and the strains are written using the Stroh formulation 
as follows: 



Qui + Ru 2 , 
R T u 1 + Tu ; 



(110) 
(111) 



The 3 x 3 m atrices Q,R and T depend on the material constants, and are defined as follows 
(jTinel . Ilflflfil ): 

Qik = Ctlfcl) Rik = Cajfc2j = Ci2k2i (H2) 

where Cijki are component s of the ela s tic st iffness tensor. Using this notation, the static equi- 
librium equations become ( Ting . 1995 . 2000l ): 



ti,i + 1 2 , 2 = Qu 11 + (R + R T )u 12 + Tu 
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(113) 



The displacement u(%i , 3 :2)1 which is a general solution of equation ()113|) . has the following form 

{5wl Il990bl : lTineLll99fih : 

u = Fg(z) + Fg(z), (114) 
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also the derivative of the displacements ui(ii,i2) and the traction t 2 (a;i,:c 2 ), can be written 
in the same form: 

t 2 (xx, x 2 ) = Lh(z) + Lh(z), (115) 

and 

u,i(xi,x 2 ) = Fh(z) + Fh(z), (116) 
where h(z) = dg/dz, F and L are constant 3x3 matrices, defined as follows: 

F = (fi,f 2 ,f 3 ), L = (l 1 ,l 2 ,l 3 ). (117) 

Note that g(z) and h(z) are analytic functions vectors with components gj(x\ + /ijX 2 ), and 
hj{x\ + /ijX 2 ) and fij are complex numbers with positive imaginary parts. According to Suol 
()1990bl ). if gj(zj) and hj(zj) are analytic functions of Zj = x\ + ^x 2 in the upper half-plane (or 
in the lower half-plane) for one fij, where fij is a complex number with positive imaginary parts, 
they are analytic for any fj,j. On the basis of this property, here and in the text that follows, the 
analysis is reduced to a single complex variable. The eigenvectors fj and the eigenvalues fij are 
simultaneously determined inserting express ion (II 141) in to equilibrium equations (|113p . so that 
they are reduced to the eigenvalue problem ( Tind . Il996l ): 



0. 



[Q + (R + R J )/i i + T^]f i 
Moreover, lj are related to fj as follows (|TineL 1 199.4 1 19961 . l20()d ): 

lj = [R T + jUjTjfj. 



(118) 



(119) 



The Hermitian surface admittance tensor, needed in general weight functions expressions ([13 



and (HU, is defined as Y = iFL -1 (jGao et allfl992h . 



Stroh form alism have b een extensively used i n analysis of interfacial cracks in anisotropic 



bimaterials by ISud (jl990bl ) and iGao et al 



(| 19921 ). Physical displacements and stress fields at 
the interface between the two materials have been derived, expressing the boundary conditions 
in terms of non-homogeneous Riemann-Hilbert problem and obtaining an algebraic eigenvalue 
involving the symmetric bimaterial matrix ()1T[) . that is solved in closed form. Recent l y, thi s 
approach has been extended and used together with Fourier transform by iMorini et al .1 (120131 ) 
for deriving singular solutions U of the elasticity problem with zero traction on the faces where 
the crack is placed along the positive semi-axis x\ > 0, and for evaluating general expressions for 
symmetric and skew-sym metric weight functions defined as traces of these functions, following 
the method illustrated in IPiccolroaz 61111 tod ). 



A. 2 Monoclinic materials 



of the stiffness tensor C. 



ijkl 



(Suo|. 


1990bl: 


Ting, 


19951 



0, employing the contracted notation 
), the three matrices Q, R and T are 
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given by: 



Cll C16 
Q = | Cie C 66 

C 55 , 




C66 C26 
C26 C 2 2 
C44 



This structure of the Stroh matrices is a consequence of the fact that for this class of anisot ropic 
materials, the redu c 
1990bl : lTineLl200d ): 



materials , the redu ced 6x6 stiffness tensor Cij possesses the following important property (jSuo 



Cl4 = Cl 5 = C 2 4 = C 2 5 = C 46 = C 56 = 0. 



(120) 



Substituting matrices Q,R and T into constitutive relations (jllOp and (jllip . we can easily 
observe that as a consequence of property (|120p . for monoclinic materials subject to two- di- 
mensional deformations in-plane strain and antiplane strain are uncoupled and also in-plane 
stress antiplane stress are uncoupled. More discuss ions and details on deco uplin g of in-plan e 
and antiplane deformations and stress are given by iHorgan fc Miller! (jl994T ) and [Ting] (|2OO0h . 
The eigenvalues problem (1118P then reduces to: 



Cn + 2ci 6 ^ + C 6 6/U 2 
C16 + (C12 + Cm)H + C 2 6V 2 





C16 + (C12 + 036^ + C26A* 2 
C 66 + 2c 2 6Ai + C22M 2 










c 55 + 2c 45 ^ + c 44 // 2 



f = 0. (121) 



Through the decoupling, the sixth order characteristic equation of this eigensystem consists in 
the product between a fourth order term correspondin g to in-plane strain and a second order 
term associated to antiplane deformations ( Tine . 19961 ): 



(cn + 2c i6 /i + C 66 /i 2 )(c 66 + 2c 26 /U + C 22 /U 2 ) - (Ci 6 + (C12 + C 66 )/i + C 26 /U 2 ) 2 = 0, 



C55 + 2c 45 ^ + C44/X = 0. 



(122) 
(123) 



Since these equations possesses real coefficient s, the roots are complex conjugates. Considering 
only eigenvalues with positive imaginary part ( Suol . 1990b ). fj,\ and \i<i are assumed to be roots 
of equation (|122p . while \iz is root of ()123p . 

The illustrated Str oh representation o f the elasticity problem is equivalent to the matrix for- 
mulation proposed by Lekhnitskii ( 19631 ). which provides alternative expressions for the eigen- 
vector matrices F and L in function of the elements of the compliance matrix S = C _1 . More 
precisely, Lekhnitskii approach gives a specially normalized eigenvector matrix F, and express- 
ing the elements of the stiffness matrix in function of the stiffness parameters it is easy to check 
that c haract e ristic equatio n derived using Lekhnitskii formulation is identical to (|122p and (|123|) 
(jSud . Il990bl : iHwul . Il993bl ). In order to obtain compact expressions for the surface admittance 
tensor Y and then for bimaterial matrices (117p and (|18p . particularly convenient for weight 
functions derivation, we assume this particular normalization for matrices F and L, reported 
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bv IHwuI (jl993bh . and we express the stiffness reduced tensor elements Qj in function of the 
elements of the reduced compliance matrix s 4 -: 



Si3S3j 
S33 



(124) 



An alternative form for characteristic equations (|122j) and (|123p is derived: 

2s' 16 // + (2s' 12 + Sg 6 )/i 2 - 2s' 26 /i + s' 22 = 0, 

/ / / n 

0. 



(125) 

S44 - 2s 45/ i + S 55^ = 0. (126) 

The hermiti an matrix Y = zFL -1 evaluated using the eigenvector normalization reported in 
Hwu ( 1993bl ) assumes the form: 



iFL 



Sn^i + M2) 

s 22 



'12 



'12 



5 22 



V 







si M2 





S11M1M2 

1 



'26 



(127) 



6 45 / 



where fi± and /i2 are roots of the equation (|125p . and [13 is solution of the equation (|126j) 
with positive imaginary part, corresponding respectively to plane str ain and antiplane strain . 
Employing the relations between fi± and and co efficients of (j!25[) (|Sud . Il990bl : iTind . Il995l ) 
and between fi^ and coefficients of fll26|) (|Sud . ll990bl ). through some manipulations the following 
form for Y is obtained: 



iFL 



I s' n Im(/ii + ix 2 ) 
*0'llPlM2 - s'12) 

V 



# H s i2 - SllMl^J 
s'uImOiijtaOZi + Pa)) 









(128) 



'44° 55 



'45, 



This general compact expression for the surface admittance tensor has been used for deriving 
explicit bimaterial matrices (|17p and (|18p and symmetric and skew-symmetric weight functions 
(fl~5l) and (fTBT) . On the basis of property (|120|) . Mode III crack propagation (antiplane shear) 
will be treated separately from Mode I and Mode II (plane strain), as for the case of interfacial 
cracks in two-dimensional isotropic bimaterials. 



References 

Arfken, G. B., & Weber, H. J. (2005). Mathematical methods for physicists, 5th edition. Elsevier 
Academic Press, San Diego. 

Barber, J. R. (2002). Elasticity. Kluwer Academic Publishers, Dordrecht. 



26 



Bigoni, D., & Capuani, D. (2002). Green's function for incremental nonlinear elasticity: shear 
bands and boundary integral formulation. J. Mech. Phys. Solids, 50, 471-500. 

Budiansky, B., & Rice, J. R. (1979). An integral equation for dynamic elastic response of an 
isolated 3-d crack. Wave Motion, 1, 187-192. 

Bueckner, H. F. (1985). Weight functions and fundamental fields for the penny-shaped and the 
half plane crack in three-space. Int. J. Solids Struct., 23, 57-93. 

Duduchava, R. (1976). Integral equations of convolution type with discontinuous presymbols, 
singular integral equations with fixed singularities and their applications to some problems 
of mechanics. Trudi Tbiliskogo Mathematicheskogo Institute/, Academii Nauk Gruzinskoi, SSR 
60, 1-135. 

Duduchava, R. (1979). Integral equations with fixed singularities. Teubner-Texte zur Mathe- 
matik. Leipzig. BSB B. G. Teubner Verlagsgesellschaft, 172 S., 10 Abb., M 16, BN 6659601. 

Gakhov, F. D., & Cherski, Y. I. (1978). Equations of the convolution type. Nauka, Moscow. 

Gao, H. (1992). Weight function method for interfacial cracks in anisotropic bimaterials. Int. 
J. Fract., 56, 139-158. 

Gao, H., Abbudi, M., & Barnett, D. M. (1992). Interfacial crack-tip field in anisotropic elastic 
solids. J. Mech. Phys. Solids, 40, 393-416. 

Gohberg, I., & Krein, M. G. (1958). Systems of integral equations on a half line with kernels 
depending on the difference of arguments. Uspekhi Mat. Nauk., 13, 3-72. 

Graciani, E., Mantic, V., & Paris, F. (2007). On the estimation of the first interpenetration 
point in the open model of interface cracks. Int. J. Fract., 143, 287-290. 

Gupta, V., Argon, A. S., & Suo, Z. (1992). Crack deflection at an interface between two 
orthotropic media. J. Appl. Mech., 59, S79-S87. 

Horgan, C. O., & Miller, K. L. (1994). Antiplane shear deformations for homogeneous and 
inhomogeneous anisotropic linearly elastic solids. J. Appl. Mech., 61, 23-29. 

Hwu, C. (1993a). Explicit solutions for collinear interface crack problems. Int. J. Solids Struct., 
30, 301-312. 

Hwu, C. (1993b). Relations among Stroh, Lekhnitskii and Muskhelishvili formulations. In Proc. 
of the 10th National Conference on Mechanical Engineering, (pp. 517-525). Oxford University 
Press. 



27 



Krein, M. G. (1958). Integral equations on an half-line with kernels depending upon the difference 
of the arguments. Uspekhi Mat. Nauk., 13, 73-120. 

Kupradze, V., Gegelia, T., Basheleishvili, M., & Burchuladze, T. (1979). Three-dimensional 
problems of the mathematical theory of elasticity and thermoelasticity . North-Holland, Ams- 
terdam. 

Lekhnitskii, S. G. (1963). Theory of elasticity of an anisotropic body. Holden-Day, San Francisco. 

Linkov, A. M., Zubkov, V. V., & Kheib, M. A. (1997). A method of solving three-dimensional 
problems of seam working and geological faults. J. Min. Sci., 33, 295-315. 

Mantic, V., & Paris, F. (2004). Relation between SIF and ERR based measures of fracture mode 
mixity in interface cracks. Int. J. Tract., 130 , 557-569. 

Mikhlin, S., & Prossdorf, S. (1980). Singular integral operators. Akademie-Verlag, Berlin. 

Morini, L., Radi, E., Movchan, A. B., & Movchan, N. V. (2013). Stroh formalism in analysis of 
skew-symmetric and symmetric weight functions for interfacial cracks. Math. Mech. Solids, 
18, 135-152. 

Mushkelishvili, N. I. (1946). Singular integral equations. Moscow. 

Mushkelishvili, N. I. (1953). Some basic problems of the mathematical theory of elasticity. 
Noordhoff. 

Pan, E. (2000). Three-dimensional Green's functions in an anisotropic piezoelectric solids. Int. 
J. Solids Struct, 37, 943-958. 

Pan, E. (2003). Three-dimensional Green's functions in an anisotropic half-space with general 
boundary conditions. J. Appl. Mech., 70, 101-110. 

Pan, E., & Yuan, F. G. (2000). Three-dimensional Green's functions in anisotropic bimaterials. 
Int. J. Solids Struct, 37, 5329-5351. 

Piccolroaz, A., & Mishuris, G. (2013). Integral identities for a semi-infinite interfacial crack in 
2D and 3D elasticity. J. Elasticity, 110, 117-140. 

Piccolroaz, A., Mishuris, G., & Movchan, A. B. (2007). Evaluation of the Lazarus-Leblond 
constants in the asymptotic model for the interfacial wavy crack. J. Mech. Phys. Solids, 55, 
1575-1600. 

Piccolroaz, A., Mishuris, G., & Movchan, A. B. (2009). Symmetric and skew symmetric weight 
functions in 2D perturbation models for semi-infinite interfacial cracks. J. Mech. Phys. Solids, 
57, 1657-1682. 



28 



Piccolroaz, A., Mishuris, G., & Movchan, A. B. (2010). Perturbation of mode III interfacial 
cracks. Int. J. Fract., 166, 41-51. 

Prossdorf, S. (1974). Some classes of singular equations. Akademie-Verlag, Berlin. 

Sinclair, J. E., He Hirth, J. P. (1975). Two dimensional elastic Green function for a cracked 
anisotropic body. J. Phys. F: Metal Phys., 5, 236-246. 

Stroh, A. N. (1962). Steady state problems in anisotropic elasticity. J. Math. Phys., 1^1, 77-103. 

Suo, Z. (1990a). Delamination specimens for orthotropic materials. J. Appl. Mech., 57, 627-634. 

Suo, Z. (1990b). Singularities, interfaces and cracks in dissimilar anisotropic media. Proc. R. 
Soc. Lond. A, 427, 331-358. 

Suo, Z., Kuo, C.-M., Barnett, D. M., & Willis, J. R. (1992). Fracture mechanics for piezoelectric 
ceramics. J. Mech. Phys. Solids, 40, 739-765. 

Ting, T. C. T. (1995). Generalized Dundurs constants for anisotropic bimaterials. Int. J. Solids 
Struct., 32, 483-500. 

Ting, T. C. T. (1996). Anisotropic elasticity: theory and applications. Oxford University Press. 

Ting, T. C. T. (2000). Anisotropic elastic materials that uncouple antiplane and inplane dis- 
placements but not antiplane and inplane stresses, and vice versa. Math. Mech. Solids, 5, 
139-156. 

Vekua, N. P. (1970). Systems of singular integral equations. Nauka, Moscow. 

Weaver, J. (1977). Three-dimensional crack analysis. Int. J. Solids Struct., 13, 321-330. 

Willis, J. R., & Movchan, A. B. (1995). Dynamic weight function for a moving crack. I. mode I 
loading. J. Mech. Phys. Solids, 43, 319-341. 

Wu, K. C. (1990). Stress intensity factor and energy release rate for interfacial cracks between 
dissimilar anisotropic materials. J. Appl. Mech., 57, 882-886. 



29 



